Dynamical symmetry breaking in transport through molecules 
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We analyze the interplay between vibrational and electronic degrees of freedom in charge trans- 
port across a molecular single-electron transistor. We focus on the wide class of molecules which 
possess quasi-degenerate vibrational eigenstates, while no degeneracy occurs for their anionic con- 
figuration. We show that the combined effect of a thermal environment and coupling to leads, 
involving tunneling events charging and discharging the molecule, leads to a dynamical symmetry 
breaking where quasi-degenerate eigenstates acquire different occupations. This imbalance gives rise 
to a characteristic asymmetry of the current versus an applied gate voltage. 

PACS numbers: 85.65.+h, 85.85.+j, 73.63.b 



Molecular electronics is a promising answer to the de- 
mand of miniaturization, reproducibility and flexibility 
of electronic devices. Starting with the pioneering work 
of Aviram and Ratner and especially with the first 
single-molecule measurement by Reed et al. @, molec- 
ular electronics has become an active research field both 

i—i 

experimentally and theoretically [3J- Still, fundamental 
questions on the peculiar nature of single-molecule junc- 
tions and their novel functionalities remain open. Among 
others they involve the unique electromechanical proper- 
ties of molecular junctions which render these distinctly 
different from, e.g., semiconductor quantum dots. 

Recent research on this particular topic has revealed a 
number of interesting effects such as shuttling instabili- 
ties 0, H, El , Franck-Condon blockade 0, @ and, more 
generally, conformational 0, Hoi Hi) or vibronic 0, 0, 
IT3. Hif signatures in the electron transport character- 
istics. In this respect, however, the role of coherences 
in transport due to quasi-degenerate levels has not been 
highlighted. 

In this Letter we consider charge transport through a 
molecular junction weakly contacted to leads. At low 
enough temperatures charging effects prevent that more 
than one excess electron at a time can populate the 
molecule. The weak coupling ensures that the poten- 
tial drop between source and drain leads is concentrated 
in the contact region, and no substantial electric field 
acts directly on the molecule. We focus on the ubiq- 
uitous, but so far poorly investigated, case of adiabatic 
vibronic potentials with more than one stable configura- 
tion as e.g. found in oligoparaphenylenes. In particular, 
we address molecules characterized by eigenstate dubletts 
in the neutral configuration and non-degenerate anionic 
states, or vice versa. This is sketched in Fig-Q] which can 
represent, e.g., the adiabatic potential energy surfaces of 
biphenyl-based molecules as a function of the dihedral an- 
gle (here playing the role of the x coordinate) between the 
phenyl rings |l5|. These molecules have been studied in 
recent transport experiments Specifically, in 

for the first time quantitatively, the connection between 
the angular conformation and the electrical conductance 
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FIG. 1: (color online) Electromechanical states of the 
molecule. Thick black lines indicate the adiabatic potentials 
Vb and Vi for the neutral (solid) and singly-charged (dashed) 
molecule. Thin green (blue) lines denote the even (odd) vibra- 
tional eigenfunctions (for the two potentials) displayed with a 
vertical shift equal to their corresponding energy eigenvalue. 
Energies are given in units of h uip and lengths in terms of 
the zero-point motion Axq = \J h/muo- The arrows mark 
examples of allowed and forbidden transitions. 



has been proved for several biphenyl-based molecules. 

We show that dynamical symmetry breaking (DSB), 
where quasi-degenerate eigenstates are differently occu- 
pied, may occur and affect transport, while in the absence 
of couplings to the leads these states are equally popu- 
lated at finite temperature, and hence the system does 
not prefer a definite parity. Temperature, on the other 
hand, should not exceed a critical value above which DSB 
is lost due to dephasing, caused by a thermal bath. Solv- 
ing the master equation for the reduced density matrix 
including coherences, we demonstrate the possibility of 
detecting DSB in the current through the molecule un- 
der different bias and gating conditions. 

To address transport in the single-electron regime, 
where only the electronic ground states |0) of the neu- 
tral and |1) of the anionic configuration are involved, we 
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describe the molecule through the Hamiltonian 



Ha 



P 

2m 



+ \0)(0\V (x) + \l)(l\V 1 (x) (1) 



in terms of the corresponding adiabatic potential energy 
surfaces V (x), Vj_(x) associated with the softest mode co- 
ordinate x. The adiabatic potentials and the related vi- 
brational eigenfunctions are sketched in Fig. ^ The neu- 
tral molecule is characterized by the double well potential 
Vq with minima located at ±A. For large enough barri- 
ers, i.e. exponentially small tunnel splitting, the low-lying 
eigenstates are organized in pairs of (quasi-)degenerate 
wave functions. The anion is modelled by a single har- 
monic potential V\ with frequency u)\, centered at the 
origin of the coordinate system (dashed line in Fig. 1). 

We express the system coupled to leads and in contact 
with a thermal bath through the Hamiltonian 



H = H s + H L + H R + H B + V T + V SB 



(2) 



The leads are described as reservoirs of non-interacting 
quasi-particles in terms of + Hr = J2k a=L rC 6 *** — 



where //l,r = ±AV/2 accounts for the effects 
of a symmetrically applied bias voltage and c\ a (c^ Q ) 
creates (destroys) an electron in the lead a. Transfer 
of electrons is mediated by the tunneling Hamiltonian 

Vt = «E feiQ =L,a(|0Xl|4 + liXOKJ . Finally, H B de- 
scribes a thermal bath of harmonic oscillators coupled 
linearly to the displacement coordinate x via the system- 
bath Hamiltonian Vsb = ^9 x ^2 q (d q + d q ), where d q and 
d q are bosonic creation and destruction operators associ- 
ated with the oscillator with energy hio q . 

Since we are interested only in the low-energy states, 
we approximate Vq(x) by two harmonic wells centered at 
±A and with frequency luq. Formally, this is achieved by 
projecting the Hamiltonian H on a smaller Hilbert space 
in which we can define the identity operator as 



l = |lXl| + |0X0|(P++7>_) 



(3) 



Here V± = J2n=o l n ' ^X 7 ^ =M are projection operators 
on the states right/left to the barrier. They are de- 
fined in terms of the shifted harmonic oscillator eigenvec- 
tors |n, ±) = exp(^fiXp/h)\n}. The condition of negligi- 
ble overlap between these eigenfunctions reads V+ V- = 
V- V+ = 0. In the Hilbert space defined by identity Q 
we obtain as effective system Hamiltonian 

Hf =!o><o| [v T {\ + S T d T )v T ) hw + 

r=± (4) 

+ |1)<1| [{\+Sd)twi + eV g -e a ] , 

where d\_ and d^ are operators for the neutral and 
charged configurations, V g is the gate voltage, and e a 



the electron affinity. Along similar lines we find for the 
effective tunneling and system-bath Hamiltonians 



vf= v £ [\o)(i\ c { a (v + +v-)+ h. c : 

fc;a=L,R 

V s c i = hgJ2(d q + dl){s\l)(l\(d + d^+ (5) 
i 

+ |0><0| J2 [V T (d T +dl + t2X)V t ]} , 



T = ± 



with g — g^2h/ (mmo) and S = y/uJi/ujQ. For transport, 
we seek stationary solutions of the Liouville equation 

&(t) = -^Tr lcads+bath {[i/ off , W(t)]} (6) 

for the reduced density matrix (RDM) a(t) := 
Tri ea ds+bath{W / (i)}, where W(t) is the total density ma- 
trix associated with the effective Hamiltonian corre- 
sponding to that in Eq. J5J. Treating the interac- 
tions V£ s and V s e g as perturbations, we rewrite Eq. © 
in terms of three contributions: a(t) = £a(t) — 
(£coh + Aim + £damp) <j{t). The coherent part of the 
differential equation has the usual form £ co h p{t) — 
— {i/h)[H§ s , <r(t)] describing the evolution of the isolated 
molecule. The coupling to the leads/bath gives rise to 
the driving/damping terms C tun a(t) and ^dampC^)- 

We now perform the following standard approxima- 
tions: we treat the leads as thermal reservoirs at equilib- 
rium temperature T, consider the coupling to the leads 
up to second order in V£ and, being interested in the 
stationary solution, neglect non-Markovian contributions 
to the equation of motion. Furthermore, neglecting off- 
diagonal elements of the RDM between anionic and neu- 
tral configurations, as well as between non-degenerate 
eigenstates, and disregarding non-energy conserving con- 
tributions, the anionic and neutral state components read 



(Aun c) 11 



(£tun C )[ 



V \2{Tf n a m V T +V T a m T^) 



a=L,R;r=± 



E Tr 

a=L,R 



r^au + crnr^J- (7) 
- 2 (I 1 hi cr oo +coor m ) V T i , 



with rate matrices for tunneling into/out of the molecule 



r fn = -7T E \m)(n,T\f(eV g +e nm )(m\n,T) 



m.n.r— ± 



E W,r){m\-T[ 



« t 



(8) 



m,n,r— ± 



Here /(e) is the Fermi function, and r Q = (2n /K)D a \v\ 2 
is the bare electronic rate, with D a the density of states in 
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Neutral state Parity 



Anion state Parity 



(a) 2<AV<4, V B =0 



(units of neo/e) 
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FIG. 2: (color online) Left hand side. Onset of the dynamical symmetry breaking with decreasing temperature (in units of 
hwo/kB) from bottom to top. The columns show the absolute value of the stationary current and the parity of the neutral and 
anionic state as a function of the gate voltage V g and the bias AV across the molecular junction, both in units of huJo/e. The 
other parameters are: 7 = 0.01 uo, ^l,r = 0.1 wo, A = 3 \J h/mujo. Right hand side. Rates scheme at the even (a) and odd (b) 
transition points (corresponding to the points marked by (a) and (b) in the upper middle panel of the figure to the left. The 
horizontal lines mark the electromechanical states of the neutral and anionic configuration with green (blue) color representing 
the even (odd) parity. Straight arrows represent transitions due to the coupling to the leads: the rate is qualitatively expressed 
by the thickness of the line. Dashed arrows show transitions lifting the blocking character of the ground states and are crucial 
for the understanding of the symmetry breaking transition. Wavy arrows represent the effect of the phononic bath. 



lead a at the Fermi energy. For convenience we assumed 
equal frequencies luq = u)\ of the neutral and anionic po- 
tential, and hence e nm = hujo(m — n); calculations for 
u>o 7^ k>i yielded no qualitative changes. Furthermore we 
shifted the electrochemical potential eV g — e a — > eV g in 
the argument of the Fermi function. The overlap matrix 
elements (m\n,r) of the vibrational states of the neutral 
and charged configurations, known as Franck-Condon 
factors, determine, together with energy resonance con- 
ditions, the transition rates © and thus the transport 
characteristics of the junction. While the Franck-Condon 
factors are fixed by the adiabatic potentials, the bias and 
gate voltages influence the resonance conditions. 

To describe relaxation and dephasing, we consider the 
case of an Ohmic bath with spectral density J(u>) = 
r«7W, where 7 is the damping coefficient ,18]. Along the 
same lines as for the tunneling term we find 



(^damp 0") j ^ 



(^damp ^) 



(ill 



-^(A> +!)[*,[,, an]], 
7muo 



(9) 



h 

"87^A^ + 0-53- V- + V- T>+) . 



where N(ui) is the Bose distribution function. Note the 
occurence of a term responsible for pure inter-well de- 
phasing and proportional to T and A 2 in the expression 
for (£ 

damp 

Since we are interested in the long-time properties, we 
look for solutions of the stationary problem £cr stat = 0. 
Given the stationary RDM cr stat we calculate the sta- 
tionary current as the trace over the system degrees 
of freedom of the left/right current operators, I stat = 
Tr roech [cr stat 7 L ]=Tr roech [a stat 7 R ] , with, e.g., 



Y, [2\o)(o\(v r Tf n + rfJv T y 



T = ± 



\l)(l\(V T T^ ut +T^Jv T ) 



(10) 



We stress that not only populations but also coherences 
of the RDM in the (±) basis contribute to 7 stat . 

In the first column of the left side of Fig. [21 we present 
the results of our calculation of the gate- and bias- 
dependent current for decreasing temperature (from bot- 
tom to top). Besides the evidence of a Franck-Condon 
blockade 0, EJ (where due to an exponential suppres- 
sion of the Franck-Condon factors (m\n, r) transport is 
blocked and the Coulomb diamonds no longer close) and 
current steps due to phononic excitations, Fig. shows 
an increasing asymmetry in the gate voltage V g with de- 
creasing temperature. To understand the origin of this 
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asymmetry, we depict in the second and third column 
of the left side of Fig. 2 the parities of the neutral and 
anionic state defined as 

N 

Pa = [(0,e,n\a\0,e, n) - (0, o,n\tr\0,o,n)] , (11) 

n=0 
oo 

Pi = f<l,2n|(7|l,2n)-<l,2n+l|(7|l,2n + l)] , 

n=Q 

where the even/odd states of the neutral molecule are 

|0,e/ O) n) = ^(|n,+)±(-l)>,-)). (12) 

In the parameter region where the current becomes asym- 
metric with respect to V g , degenerate states are differ- 
ently populated, i.e., a dynamical symmetry breaking 
(DSB) in the occupation of the even/odd states occurs. 
In the regions of defined parity (warm colors for even, 
cold for odd) the system exhibits spatial coherences in 
the (±) basis, cf. Eqs. In contrast, we find that in 

the even/odd basis coherences are decoupled from popu- 
lations and vanish for the stationary density matrix. The 
onset of this even/odd parity regions in the gate-bias volt- 
age plane can be understood in terms of transition rates 
between the different vibrational states of the molecule 
(see the right side of Fig. [2I> - At zero bias, energy conser- 
vation prevents any transition from the neutral (even or 
odd) or anionic ground states to occur. These states are 
thus blocking states. In particular for positive gate volt- 
ages the molecule is neutral and dephasing ensures equi- 
libration of the even and odd populations of the degener- 
ate ground state (P — 0) . At negative gate voltages the 
molecule is charged and dissipation ensures relaxation 
predominantly to the anionic ground state (Pi « 1). 

The situation changes for combinations of bias and 
gate voltages that allow current to flow: the electrome- 
chanical system is maintained out of equilibrium by the 
applied voltage. Now the parity of the anionic state fluc- 
tuates due to population of the higher excited states, and 
the distribution of the populations no longer corresponds 
to the thermal distribution in equilibrium with the bath. 

The sign of the neutral-state parity depends on the par- 
ticular rate configuration. Though a quantitative predic- 
tion relies on the solution for er stat , it is possible to under- 
stand the sign of the neutral state parity with simple ar- 
guments. In the right side of Fig. [21a representation of the 
electromechanical rates for two specific cases is reported. 
Despite the complexity of the scheme, only few lines are 
crucial: these are the two red dashed arrows representing 
the rates which lift the blocking character of the neutral 
ground states. In the case (a) the neutral odd ground 
state can be depopulated resulting in an even parity of 
the stationary state. In case (b) the neutral even ground 
state is involved leading to an odd parity of the stationary 
state. These unblocking rates represent the bottle-neck 



in the DSB and are competing with the dephasing gen- 
erated by the heat bath (wavy line joining the even and 
odd sector of the neutral state) . The intensity of the un- 
blocking rates depends on the Franck-Condon factors of 
the vibrational wave functions involved. A comparison 
with the dephasing rate allows to estimate the transition 
temperature T tI as fcsTtr = fi^or u b/(4A 2 7) , where F u b 
is the unblocking rate (i.e. the depopulation rate of the 
relevant neutral ground state). The even transition (a) 
involves a larger r u t,: it is more robust and occurs at 
higher T (compare first and second panel in the central 
column of the left side of Fig. |5J ■ 



To summarize, we analyzed the dynamics of a molecu- 
lar junction in the single-electron transport regime. The 
molecule possesses quasi-degenerate vibrational eigen- 
states in the neutral configuration and no degeneracy 
in the anionic case. Tunneling processes charging and 
discharging the molecule preserve the parity of the wave- 
functions. As a consequence, unequal occupation of de- 
generate molecular neutral states occurs. An explanation 
of this effect in terms of unblocking rates is given. 
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M. Cfzek for inspiring discussions. 
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